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ABSTRACT 

Three Galactic O-type stars belong to the rare class of Of?p objects: HD 108, HD 191612, and 
HD 148937. The first two stars show a wealth of phenomena, including magnetic fields and strong 
X-ray emission, light variability, and dramatic periodic spectral variability. We present here the 
first detailed optical and X-ray study of the third Galactic Of?p star, HD 148937. Spectroscopic 
monitoring has revealed low-level variability in the Balmer and Hen A 4686 lines, but constancy 
at He I and C in A 4650. The Ha line exhibits profile variations at a possible periodicity of ~7d. 
Model atmosphere fits yield T eff = 41000 ± 2000K, \og(g) = 4.0 ± 0.1, M sph < lO^Moyr" 1 
and a surabondance of nitrogen by a factor of four. At X-ray wavelengths, HD 148937 resembles 
HD 108 and HD 191612 in having a thermal spectrum dominated by a relatively cool component 
(fcT=0.2keV), broad lines (>1700kms _1 ), and an order-of-magnitude overluminosity compared 
to normal O stars (log[L^ nabs /L B oL] 6). 

Subject headings: stars: individual: HD 148937 - stars: early-type - X-rays: stars - X-rays: individual: 
HD 148937 
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Introduction 



HD 191612 to be a bin ary, though with a d iffer- 



The Of?p category was defined by IWalborn 



(|1972l ) to gather peculiar stars displaying strong 
Cm emission lines around 4650 A, i.e. of an in- 
tensity comparable to that of the neighbouring 
Nni lines. Three stars of our Galaxy belong to 
this class: HD 108, HD 191612, and HD 148937. A 
few others ha ve now been identified in the Magel- 
lanic Clouds (|Hevdari-Malaveri fc Melnick 1992; 
IWalborn et al. Il2000t iMassev fc Duffy Il2001f) . 

In recent years, these peculiar objects have at- 
tracted quite a lot of attention. In 2001, it was 
discovered that HD 108 presents dramatic line- 
profil e variations of its Cm, He I and Hydrogen 



lines (N aze et al. Il200lh . These variations seemed 



recurrent with a timescale of app r oximately 50-60 



years (see also lNaze et al. 1120061 ) . IWalborn et al 



(|2003l ) ;hen reported a similar phenomenon for 
the spectrum of HD 191612. However, the period 
was much shorter (only 538d) an d appeared cor 



relate d with photometric changes (jWalborn et al 



20041 ). A magneti c field was subseque ntly iden- 



tified for this star ( Donati et al. 20061) . with the 



538 d proposed to be linked to the rotation period. 
In this scenario, HD 191612 would be an oblique 
magnetic rotator, a somewha t evolved version of 



9 1 Ori C (jDonati et al. Il2006h . The measured pe- 



riod would correspond to the rotational period of 
the star, which has been slowed down by the mag- 
netic field, and larger emissions would be detected 
when the magnetically-confined disk is seen face- 
on while there are reduced or no contamination by 
emission when the disk is seen edge-on. 

However, the analysis of a dedicated XMM- 
Newton observing campaign did not reveal the 
dominant signature of a magnetically-confined X- 
ray emitting plasma: apart from a clear over- 
luminosity, both HD108 and HD 191612 dis- 
play characteristics typical o f massive O-type 



stars (broad lines, low kT; see iNaze et al. 112004 . 
20071) . Optical observations have further shown 
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ent period (1542d, see iHowarth et al. 1 120071 ). A 



magnetically-confined disk therefore appears in- 
sufficient to explain the full behaviour of those 
Of?p stars; an additional source of relatively soft 
X-rays is needed but its origin is still unknown 
(is it somehow related to binarity, like e.g. col- 
liding wind emission, or is it a more exotic phe- 
nomenon?) 

In comparison, only little is known about 
HD 148937, and no thorough study of the star has 
been undertaken until now. The most interest- 
ing information comes from its environment: the 
star is surrounded by a bipolar nebula, NGC 6164- 
6165, that displays a very symmetrical geometry. 
Due to its anomalous chemical abundances, this 
nebula is thought to have formed through an erup- 
tion of the star, maybe simil ar to that of 77 Carina e 



in the nineteenth century ( Dufour et al. 19881 ). 



making HD 148937 a candidate Luminous Blue 
Variable object. 

The aim of our study is to assess the properties 
of HD 148937, in the visible and X-ray domains, 
and to compare it to the other two galactic Of?p 
stars. This paper is organized as follows. Section 2 
describes the observations while Sections 3, 4 and 
5 present the results of the available photomet- 
ric datasets, the dedicated optical spectroscopic 
campaign and the XMM-Newton observations, re- 
spectively. Section 6 gives our conclusions. 

2. Observations and Data Reduction 
2.1. Visible domain 

Low-resolution spectroscopy of HD 148937 was 
gathered between 2003 and 2006 via SMARTS 
(Small and Moderate Aperture Research Tele- 
scope System) at the CTIO 1.5m telescope 
equipped with the RC spectrograph. For each 
observation, three exposures of one to two min- 
utes were combined, and the final signal-to-noisc 
ratio (S/N) in the continuum was generally about 
100. Three different gratings and/or settings were 
used: one covering the range 5650-6800 A with a 
resolution of 2. 8 A the second 4000-4900 A with 
a resolution of 2.0A, and the last 4065-4700 A 
with a resolution of 1.6A. All the reductions were 
performed using the IRAF and MIDAS software. 
The normalization was done by fitting polyno- 
mials through carefully chosen continuum win- 
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dows. The observations were relatively evenly 
distributed between the months of February and 
October. During the 4 years of the observing 
campaign, at least one bluc+rcd spectrogram pair 
was taken each month of this visibility window. 
In total, we collected 53 red spectra and 56 blue 
ones. 

Additional high-resolution spectra covering the 
whole optical range (3750-9200 A) were collected 
using the Fiber-Fed Extended Range Optical 
Spectrograph (FEROS), an echelle instrument 
mounted at the MPG/ESO 2.2-m telescope at 
La Silla. The resolution of FEROS is 0.1 A at 
4800 A (thus R=48000); the exposure time was 
15 min, leading to a signal-to-noise ratio of mini- 
mum 300. The spectra were all taken during the 
month of May (2003 May 24 and 25, 2004 May 6 
and 2005 May 19). These data were reduced by 
the observers, E. Gosset and H. Sana, using an im- 
prove d version of the FEROS context ( Sana et al. I 



2003) working under the MIDAS environment. Af- 



ter correcting for the blaze using flat fields, the 
different orders were normalized individually us- 
ing polynomials of order 4-6. 

Other high-resolution observations can be 
found in the archives. A first one was taken 
with FEROS on 2005 June 26 and was kindly 
made available to us by the PI, J.-C. Bouret 
(S/N=200). Another one, taken on 1995 April 
15 with University College London Echelle Spec- 
trograph (UCLES), was found in the AAT archive. 
It consists of three exposures of 3-5min with S/N 
of about 140 that provide a spectrum covering 
the wavelength range from 3800 to 7400 A with 
a resolving power of 36000. The reduction of this 
observation was done in a classical way using the 
IRAF and MIDAS software. Finally, the star was 
also observed on 2002 February 27 with the Ultra- 
violet and Visual Echelle Spectrograph (UVES), 
in the framework of the UVES POP programme. 
The resolution was about 0.05A, and the signal- 
to-noise ratio ^250. The reduced and merged 
UVES spectrum was downloaded from the pub- 
lic ESO databas^E Only a few selected regions 
were extracted from the UVES spectrum. All 
these spectra were corrected for the blaze using 
flat fields and then normalized using low-order 
polynomials. 



2 http: / / www.sc.eso.org/santiago / uvespop 



2.2. X-ray range 

On 2001 Feb. 25 (Rev. 0223, PI R. Smith), 
HD 148937 was observed twice consecutively with 
XMM-Newton, for a total exposure time of 30ks. 
The star is associated with sources J163352.4- 
480640 and J163352.3-480641 in the 1XMM cat- 
alog. We retrieved the two datasets from the 
XMM-Newton public archives, in order to per- 
form a thorough analysis. For these observa- 
tions, the three European Photon Imaging Cam- 
eras (EPICs) were operated in the standard, full- 
frame mode, except for the second pn dataset 
which was taken in the extended full frame mode. 
A medium filter was used to reject optical light. 
We processed these archival data with the Science 
Analysis System (SAS) software, version 7.0. Af- 
ter the pipeline chains (tasks EMPROC, EPPROC 
and RGSPROC), the data were filtered as recom- 
mended by the SAS team: for the EPIC MOS 
(Metal Oxide Semi-conductor) detectors, we kept 
single, double, triple and quadruple events (i.e. 
pattern between and 12) that pass through the 
#XMMEA.EM filter; for the EPIC pn detector, 
only single and double events (i.e. pattern be- 
tween and 4) with flag=0 were considered. To 
check for contamination by low-energy protons, 
we have further examined the light curve at high 
energies (Pulse Invariant channel number> 10000, 
E>10keV, and with pattern=0). No single, indi- 
vidual flare was detected during the observations 
but the background level was rather high and quite 
variable during the whole exposure. Further anal- 
ysis was performed using the FTOOLS tasks and 
the XSPEC software v 11.2.0. 

3. Photometry 

HD 148937 was included in the "New Catalogue 
of Suspected Variable Stars", under the entry 
NSV7808. In this catalogue, it presents a V mag- 
nitude between 6.71 and 6.81 mag. However, the 
variabili t y stat us of this star is still under debate. 
Balona ] (|l992h found a possible dimming of the 
star b y about 0.01 mag over a few w eeks time. Fol- 
lowing |ro^Gendere^fiE5ZI ( 1989h . the star rather 
presents a constant luminosity in V , although with 
a dispersion of 0.005 mag in April 1983, and pos- 
sible small color variations (^0.002 mag). Com- 
paring with older data, the same authors also 
noted that HD 148937 might have been bluer and 
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brighter in the l ate ei ghties than in 1960. Finally, 



van Genderen I (|200ll ) considers HD 148937 as a 



possible candidate SDor variable, but with only 
few indications for its SDor status. 

The three galactic Of?p stars were observed 
by the Hipparcos satellite^. We have downloaded 
the individual photometry measurements from the 
Hipparcos database, including the Tycho data in 
the B and V filters and the Hipparcos broad band 
filter Hp. After discarding the measurements 
with non-zero flags to get rid of possibly prob- 
lematic data, we have analyzed the photometry 
of HD 148937 (=HIP 81100, see Fig. QJ. The star 
displays a H p magnitude varying between 6.79 
and 6.83 mag, but the data are relatively sparse. 
Nevertheless, we have determined the reduced \ 2 
when fitting the data by a constant luminosity, cal- 
culated the autocorrelation for each dataset, and 
performed a period search. No significant change 
or periodicity was detected: the variable/S Dor 
status of HD 148937 still awaits confirmation. 

4. Visible spectroscopy 

The full visible spectrum of HD 148937 is shown 
in Fig. [5] Compared to recent, quiescent spec- 
tra of HD 108, that of HD 148937 presents ap- 
preciably weaker absorption lines of NmAA4510- 
4534 and Hen. The faint Sim, On, and Cn 
emission lines seen in HD 108 do not exist here, 
OnA4705 being rather in absorption. The two 
SilvAA4088,4116 lines are however clearly in emis- 
sion, whereas they were in absorption in the spec- 
tra of HD 108. We may also note the important 
strength of the HeiA6678 emission line, compared 
to Ha, in the spectrum of HD 148937. In addition, 
no P Cygni profile is seen for the He iAA4471,4713 
lines nor for the H7 or H/3 lines (all are pure ab- 
sorption). However, their profiles are not symmet- 
ric, with the red wing being much steeper than 
the blue one, indicating a possible contamination 
by some emission, most probably coming from the 
wind. 

The spectral type of HD 148937 appears rather 
early. In the high-resolution spectra, a visual com- 



3 The database of the All-Sky Automated Survey (ASAS) 
also contains some data on HD 148937, but they are not 
usable: as is often the case for relatively bright stars, even 
the photometry flagged as good appears somewhat erratic, 
probably because of saturation problems. 



parison of the He IA4471 with He IIA4542 lines and 
of HeiA4026 with HeiiA 4200 favors a 06 spec- 
tral t ype (see Fig. [2] and Walborn fc Fitzpatrickl 
2000). From the measurements of the equivalent 
widths (EW) of the former pair (see below), the 
high-resolution spectra again favors 06, whereas 
the low-resolution SMARTS data rather give 05.5, 
but close to the limit between the 05.5 and 06 
types (see Fig. [3]). The star has been previ 



ously classified as 06 (IMac Connel l & Bidclman 



Humphreys 111975^ . These variations, if real, may 
indicate a similar behaviour as seen in HD 108 and 
HD 191612. However, these older spectral type 
determinations were reported more or less at the 
same epoch, and the differences could rather be 
due to different interpretations of the same spec- 
trum. 

4.1. Radial velocities and equivalent widths 

We investigated our 4 years of low-dispersion 
spectra in order to search for large, monthly to 
yearly variations. The radial velocity (RV) of the 
main spectral lines was measured by fitting a gaus- 
sian to the top (resp. bottom) of the emission 
(resp. absorption) lines. The EW was determined 
by integrating the profile in a window of about 
8A centered on the rest wavelength (such a large 
window was necessary because of the low reso- 
lution of the SMARTS data). The mean of the 
RVs and EWs and dispersion around this mean 
are presented in Table [T] Note that the RVs of 
He IA4471 appear systematically bluer than those 
of He nA4542 because of the contamination of the 
red wing of the Hei line. 

The RV and EW of the diffuse interstellar bands 
(DIBs) near Ha and Civ A 5812 were also mea- 
sured. To get rid of possible remaining calibration 
problems, the observed RVs of these strong and 
rather narrow features were used to correct the 
RVs of the main lines, assuming the rest wave- 



19761 [Penny et al. 1 119961 : iGarrisson et al. I Il977l) 
Q6.5 dWalborn Ill972h and Q7 (IConti et al. Ill977 : 



lengths from iHerbigl (|1995l ) and that the average 
DIB RVs of the FEROS data were the actual ones 
(-lS.ekms^ 1 for the DIB at 5780A, -Sikms" 1 
for the DIB at 6614A). However, only a slight im- 
provement is seen in the dispersion of the RVs after 
this correction. 

The dispersion of the RVs appears rather con- 
stant among the measured lines (about 10km s _1 ). 
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The mean RVs from the high-resolution data 
(FEROS, UVES, UCLES) are compatible with 
the SMARTS values within the error bars (see Ta- 
ble Q] - note there is no dispersion for UVES and 
UCLES data since there is only one exposure). 
It is however worth noting that HeiA4471 and 
Ha, two lines affected by wind emission, present 
a higher dispersion of the RVs in the high-quality 
FEROS datasets. Our measurements are also in 
agreement with the RVs found in the literature 
(IConti et al Ill977t lAugensen Ill985l) . In addition, 



Conti et al. I (|1977l ) mention ed the compa t ibilit; 



of their results with those of lAbt fc Biggs I (1972) 
and their average value of the velocit y is also sim- 



i lar, within the errors, to those o f IWesterlund 
(|l96ll) ; lBuscombe fc Morris! (|l960h . 



From spectroscopic data, we have thus estab- 
lished a dispersion of lOkms" 1 (=ctrv) on the 
radial velocities of HD 148937. Furthermore, no 
clear signature of binarity, like e.g. the presence 
of double-line profile, was detected in our data. 
Therefore, we conclude that no indication of bi- 
narity could be found on both the long-term (a 
few years) and short-term (a few days) ranges 
probed by the observations. This result agrees 
with the conclusion reached bv lConti et al. I ( 1977 ) 
and Garmanv et al. ( 1980h on the constancy of 
the RVfl 

However, it should be noted that the impli- 
cations of the lack of significant RV variations 
in terms of the multiplicity of the star depend 
strongly on what assumptions we make about 
the orbital parameters of a putative binary. In- 
deed, unfavourable values of the orbital incli- 
nation, the mass ratio and the orbital period 
or a combination of several of these parameters 
could lead to rather small radial velocity varia- 
tions. In this context, it is worth noting that RV 
variations similar to those found for H D 191612 
(jNaze et al. 1120071 : iHowarth et al. 112007b would be 
below o ur detection threshold. The simple ap- 
proach of lGarmanv et alTI (|l980h enable us to eval- 
uate the probability that we have missed a spec- 
troscopic binary as a result of a low orbital inclina- 
tion and for a given period and mass ratio. Taking 



Using the same data as these authors, Hutchings (1976) 
had suggested the star to be a binary of period 9d, but this 
was never confirmed later. However, a probl em was iden- 
tified with some of the photographic plates (Conti et al. 
I1977T I. 



the rather conservative assumption that the semi- 
amplitude of the radial velocity curve should be 
smaller than 2 x ctrv, and assuming a circular or- 
bit and a mass of 55 M Q for the Of?p star, we 
find that there is a less than 1.7% probability that 
HD 148937 is a binary with a period of 100 days 
or less and a mass ratio of 1. If we allow a mass 
ratio of 5 (i.e. a companion of spectral type ^B2, 
the probability increases to 6 and 24% for peri- 
ods shorter than 15 and 100 days respectively. Of 
course, the probabilities of missing a binary in- 
crease drastically if the orbit is no longer assumed 
to be circular and this would be especially relevant 
for a long-period binary system. 

Looking at the EWs, three lines, He n A 4686, 
H/3, and Ha, present clearly deviant dispersions, 
suggesting intrinsic variability. Fig. |4] compares 
the RVs and EWs of Ha and the neighbouring 
DIB, presented with similar scales: the larger dis- 
persion of the Ha EWs is obvious. We might note 
that these lines are the most sensitive ones to the 
wind. 

4.2. Line profiles 

To investigate further the variability found pre- 
viously, we performed another tes t using the Tem- 
poral Variance Spectrum (TVS, iFullerton et al. 
19961 ). No large, significant variation was found, 



except for the aforementioned variable lines and 
He I A 6678. Values close to the significance thresh- 
old are however also found for N III AA 4634,4641 
(see Fig. EJ. 

Fig. [6] displays the line profile variations of Ha 
in the SMARTS data. Note that the changes occur 
on relatively short timescales, as underlined by our 
high-resolution data. In fact, two FEROS spec- 
tra taken on two consecutive nights do show the 
small scale variations in the He II A 4686 and Hy- 
drogen Balmer lines (see right side of the same fig- 
ure for H/3) . Similar changes of the same lines are 
also found when comparing with the other high- 
resolution data. The variability is confirmed by 
looking at published tracings of the spectrum of 
HD 1489 37: in a spectrum taken in 1991 and re- 
ported in Not a et al. ( 19961 ). the H/3 and H7 lines 
seemed to consist of two absorption components 
of equal strength in 1991, whereas the blue one 
is much stronger in our 2002-2003 high-resolution 
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spectra|f|. 

A period search was made on the RV and EW 



measu rements using the techniques of Heck et al- 



ii 1985, see remarks of 



Gosset et al. I [2001) and 



Lafler fc Kinman I (|1965l ). No significant period 



was detected, except maybe for Ha, where a peak 
at about 7d seems to slightly stand out above the 
noise. We also performed a 2D Fourier test on 
the SMARTS spectra themselves, in regions cen- 
tered on the Ha, H/3, and He n A 4686 lines. The 
periodogram was then averaged on the wavelength 
interval (see the result in Fig. [7J. This time, with- 
out ambiguity, a period of 7.031±0.003 is clearly 
detected for Ha. When folded with this period, 
the RVs and EWs of Ha actually follow a regular 
pattern (Fig. [8]). This frequency is also present in 
the periodograms of H/3 and He n A 4686, though 
with a reduced intensity. A short-term moni- 
toring, preferentially done with a high resolution 
spectrograph, should better constrain the proper- 
ties of this periodic phenomenon. 

On the contrary, other lines, among which He I 
and C in A 4650, appear remarkably constant for 
an Of?p star. In that regard, the 1995 UCLES 
data could be a carbon copy of the latest FEROS 
data. Even older observations, taken by J.-M. 
Vreux in June 1974 and B. Westerlund in July 
and August 1974, show similar features, with- 
out any striking differences. The only reference 
to a possible change o f these lines can be found 



in IWesterlund I (|1961I ): 'Practically all spectra 
have He II A 4686 in emission, most of them also 
N in A 4641 and a few also Cm (?) A 4651' - with- 
out having access to the actual plates, it is diffi- 
cult to judge what happened that year, but maybe 
the crude photographic plates might be to blame. 
Without further information, it is difficult to judge 
the significance of this report; all we can say is that 
such a behavior is not seen in intensive, higher 
quality, subsequent data. 

4.3. Physical parameters 

The visible spectrum can also be used to 
derive the physical parameters of HD 148937. 
For example, the line width reflects the pro- 



5 We might also note that from objective prism data, 
iMac Connell fc Bidelman I (Il976f ) reported that the H/3 line 
of HD 148937 was 'filled in' but without a tracing, we can't 
compare their result with ours in detail. 



jectcd rotation velocity of the star. Avoiding 
the variable Balmer Hydrogen and Hen lines 
and the contaminated He I lines, we applied 
the F ourier method (see ISimon-Diaz fe Herrero 



2007, and references therein) on the metal lines 
of ClV A 5812 and OmA5592. This results in 
usin(z) of SSkms" 1 for the former and 45kms~ 1 
for the latter, and we therefore conclude that 
wsin(z) > 45kms _1 . This value is much lower 
than those found in the literatur e (200km s -1 , 

" 19961: 



Conti fc Ebbets 1977; 92km s _1 , iPennv 



76km s , iHowarth et al. I [1997]) , but this is not 



surprising since the other methods do not permit 
one to disentangle the rotational broadening from 
other broadening mechanisms ( i.e. macroturbu- 



lence, ISimon-Diaz fc Herrero 1120071 ) 



An average of our FEROS spectra was fitted 
using CMFGEN (|Hillier fc Milled 1 19981 for a de- 
scri ption of the models and of the method, see 
also lMartins et alj|2005bt iBouret et al J 120051) . as 
shown in Fig. [9j Note that, if the central part of 
the Balmer lines is contaminated by emission, the 
wings remain sufficiently clean to be used for grav- 
ity estimates. The best fit gives T e // = 41000 ± 
2000K, log( 5 ) = 4.0 ± 0.1, R/Rq = 15.0 ± 2.5 and 
log(Z/L Q ) = 5.75 ± 0.1 (assuming He/H of 0.08 
and a distance of 1.38 kpc, see below). Such val- 
ues are comparable t o those of 05-6 V/III stars 



([Martins et al. 1120051 ) 



The abundance of nitrogen appears to be 
N/H=3 x 10" 4 (with an uncertainty of 40%), 
which for a s olar reference value of 7xl0~ 5 
(lAsplund 1120051) corresponds to an overabundance 
of a factor 4 compared to the Sun. This overabun- 
dance indicates that the star is already chemically 
evolved, showing products of the CNO cycle at its 
surface. This result is compatible with abundance 
estimates in the surrounding nebula, thought to 
consist at least partly of material ejected by the 



star (0.7 dex enrichment in N, see iDufour et al. 
19881) . 



In addition, we estimated the mass-loss rate 
from the P Cygni profiles of an archival IUE spec- 
trum. A rather high clumping factor (/=0.01) was 
required to prevent N iv A 1720 (and to a lesser ex- 
tent O v A 1371) fro m being too st r ong. T his value 
was also found by IBouret et al.1 (2005) in their 
study of two 04 stars. With this clumping fac- 
tor, N iv A 1720 and Si iv AA 1393,1403 are reason- 
ably reproduced for a mass loss rate lower than 1— 
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2xl(T 7 M Q yr- 1 . However, Nv AA 1238,1242 and 
C iv AA 1548,1550 remain too strong as long as M 
is larger than ~ lO _8 M0yr _1 . Given the evidence 
for non spherical emission in the star (see below), 
we refrained from going into too much detail in 
this analysis. One can simply conclude that a 
mass-loss rate of 10 _7 MQyr -1 (corresponding to 
M unc i = 10 _6 M Q yr _1 ) is a conservative upper 
limit on the mass loss rate of HD 148937. Note 
that the prefered value for the terminal velocity 
amounts to 2600km s -1 . 

The fit appears rather good (see Fig. [9|), ex- 
cept for two caveats. First, we note that all the 
observed He I and He n lines appears to be stronger 
than in a model with He/H=0.1, but reducing this 
ratio to 0.06 does not completely solve the prob- 
lem. A similar difficulty was unc overed when fit- 



ting the spectrum of HD 191612 (jWalborn et al. 



2003) and the origin of this discrepancy remains 
unclear at the moment (dilution by an hidden 
companion? contamination by emission?). In- 
creasing the slope of the velocity field (the clas- 
sical (3 parameter) does not help to weaken sys- 
tematically all He lines. Second, it is impossible to 
reproduce the visible emission lines with the mass- 
loss rate and terminal velocity estimated from the 
UV P-Cygni profiles. The Balmer and He n emis- 
sions are narrow and cannot be explained (only) 
by a spherical wind emission as derived in the UV. 
Therefore, the wind should consist of two compo- 
nents: a spherical one (the only one CMFGEN is 
able to reproduce) and a non-spherical one where 
the narrow Balmer emissions arise (e.g. a disk). 

4.4. Summary 

Contrary to HD 108 and HD 191612, HD 148937 
does not display spectacular line changes for H, 
He I and C in A 4650. However, it does show small- 
scale variations of the He II A 4686 and Hydrogen 
Balmer lines. Looking closely at the Ha line, a pe- 
riod of 7.031d is detected on our low-resolution ob- 
servations, but requires confirmation as the tem- 
poral sampling was not optimized for such short 
timescales. 

The two other Galactic Of?p stars actually 
also display some small-scale changes, in addi- 
tion to the main, large variations. For HD 108, 
it can be easily spotted when comparing obser- 
vations taken during one obs erving run ( see e.g. 
Fig. EH reproduced from Naze 120041 ). For 



HD 191612, iHowarth et al. I (|2007l) found 'small- 
amplitude variability ... on timescales longer than 
a few days'. It is not yet known if these changes are 
stochastic or periodic, as the observations were not 
optimized to detect such short-term variations. 

Although of small amplitude, the phenomenon 
observed for HD 148937 could still be related to 
those observed in the other Of?p stars. On the one 
hand, the variations of HD 148937 could represent 
a scaled-down version of those observed for HD 108 
and HD 191612, since these changes are similar 
in character. As the H and Hen A 4686 varia- 
tions were clearly the largest for these objects, it 
is possible that the He I and C in changes would 
then remain undetected for HD 148937. The lower 
amplitude of the variability for HD 148937 might 
be linked to a relaxation of the system following 
the eruption that gave rise to the bipolar nebula 
NGC6164-5, or could be related to a lower an- 
gle between the rotation and magnetic axis (in 
the case of the magnetic oblique rotator model, 
as proposed for HD 191612^). On the other hand, 
it is equally possible that these small-scale changes 
are similar to the short-term variations seen in 
HD 108 and HD 191612. It would then remain to 
explain their periodicity in HD 148937, and see if 
the small-scale changes are periodic for the other 
two objects. In this case, larger variations could 
be present in HD 148937, but only with a much 
longer timescale than for HD 108 since we see no 
significant changes of Hei and Cm between 1974, 
1995 and 2002-5 data. 

5. X-ray data 

The X-ray emission of HD 148937 was first dis- 
covered with the Einstein satellite. The star ap- 
pears as source 2E1630. 1—4800 in the Einstein 2E 
catalog (Harris et al. 1994). It has an IPC count 



rate of 0.100±0.004 cts s" 1 (lHarris et al. 
or 0.114±0.004 cts s" 1 dChlebowski et al. 



1994) 
1989h . 



More recently, it was observed by the ROSAT 
satellite, notably during the All-Sky Survey, where 
it appears as source 1RXS J163352. 2— 480643 (also 
known as 1RXP J163352.9-480635), with a PSPC 



6 In this case, the 7d period would reflect the rotation period 
of the star. The magnetic field should therefore be weaker 
than for HD 191612 since, although the star is old enough 
to have undergone eruptions, the field has not slowed the 
rotation to very long timescales as found for HD 191612. 
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count rate of 0.23±0.03 cts s _1 . ROSAT has also 
observed HD 148937 during three pointed observa- 
tions (collected under the id # RP900379). The 
first observation was made in September 1992 dur- 
ing 3.35 ks, the second was taken in March 1993 
during 10.14 ks and the last one in August 1993 
for 5.36 ks. After downloading these data from the 
archives, we estimated PSPC count rates for these 
observations to 0.197±0.008 cts s _1 , 0.189±0.004 
cts s -1 , and 0.203±0.006 cts s -1 , respectively. 

More recently, the star was observed by XMM- 
Newton, whose high sensitivity permits a deeper 
analysis of its high-energy emission. The EPIC- 
MOS spectra of HD 148937 were extracted over a 
circular region with radius 50"centred on the star; 
the background was extracted in the surrounding 
annular area of outer radius 75". For EPIC-pn, 
the radius of the source region was limited to 37'.'5 
in order to avoid a nearby gap; we then use as 
background a nearby circle devoid of sources. The 
EPIC spectra, binned to get a minimum of 10 cts 
per bin (i.e. S/N>3 in each bin), are shown in Fig. 
1111 where pn data is presented in green, MOS1 in 
black and MOS2 in red. 

The spectrum appears thermal, and the Fexxv 
line at 6.7keV is present. We have adopted a dis- 
tance of 1.38 kpc (based on the membership i n the 
Ara OBla association, see iHumphrevd 19781 ) and 
an interstellar column de nsity Nh of 4x 10 21 cm~ 2 
(jDiplas fc Savage 1994). [n the spectral mod- 



elling, we do not allow the absorption to go below 
this threshold. In addition, since both datasets 
gave consistent results, within the errors, we fi- 
nally fit all data simultaneously. The results are 
very similar to those of HD 108 and HD 191612: 
fits using only one thermal component were un- 
acceptable (x 2 > 2) whereas the sum of two ab- 
sorbecQ, optically thin e quilibrium plasma models 
(mekal, iKaastral Il992l ) is much better. In the 
latter case, two solutions with similar residuals 
are found: one with temperatures kT around 0.6 
and 2 keV and negligible absorbing columns; the 
other with lower temperatures (0.2 and 1 keV) and 
larger columns (iV H ~ 0.5 and 1 x 10 22 cm -2 ). 
The second is slightly better and also favored 
by the results of a fit by a differential emission- 



7 Since these two thermal components could arise in differ- 
ent regions of the wind, we allow them to be absorbed by 
independent column densities (-/V H ). 



measure model (DEM, cGpmekl, iLemen et al. 
19891: ISingh et al. I fl996h . The low-temperature 



peak is generally explained by instrinsic wind- 
wind shocks, whereas the high-temperature fea- 
ture could be related to colliding winds in a bi- 
nary or to shocks in a m agnetically confined wind 



(Zhe kov fc Palla 112007 ). If we assume the latter 



origin, it is important to note that the importance 
of this second peak is much lower than for well- 
known magnetic system such as 0Ori 1 C or rSco 
(jZhekov fc Palla 1120071) . 

As HD 148937 is brighter and closer than the 
two other galactic Of?p stars, the noise on its X- 
ray spectrum is reduced. This helps to spot a small 
deficiency of the aforementioned model at low en- 
ergies. The fit can be improved by the addition 
of a third mekal component. Two solutions are 
again found, each corresponding to the two solu- 
tions mentioned above plus a third cooler compo- 
nent. Although the 'hot' solution (kT of 0.2, 0.6 
and 2keV) presents a better x 2 than the 'cool' one 
(kT of 0.1, 0.2 and IkeV), it appears worse at low 
energies. 

The results of these fits are reported in Table 
where the unabsorbed fluxes /x nabs (in the 0.4- 
lO.keV range) are corrected only for the interstel- 
lar absorbing column. For each parameter, the 
lower and upper limits of the 90% confidence in- 
terval (derived from the ERROR command under 
XSPEC) are noted as indices and exponents, re- 
spectively. The normalisation factors are defined 
7 / n e nndV , where D, n e and hh are re- 



us 



10" 



irrD 



spectively the distance to the source, the electron 
and proton density of the emitting plasma. 

U sing typical color s and b olometric corrections 
from Martins fc Plez 1 (|2006l) and the optical prop- 
erties of HD 148937 (V=6 .728, B - y=0.343, 



Maiz-Apellaniz et al. I 120041 ) . the bolometric lu 



minosities amounts to Tbol ~2x 10 3y erg s , a 
value in agreement with the results of model atmo- 
sphere fits. The X-ray luminosities L™ abs , evalu- 
ated from the EPIC data in the 0.5-10.0 keV range 
and corrected for interstellar absorption, is ~ 2 x 

10 33 erg s -1 , resulting in log(L™ ab 7T B OL) 6. 

Again, as for HD 108 and HD 191612, the value 
of this ratio is much larger (about 8 times here) 



than that of the 'canonical' relation (jSana et al. 
20061 ). In this context, it is interesting to note 



that for the 2T fit, only 30% of the unabsorbed 
flux comes from the high-temperature component 
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(in comparison, it is 25% for HD 108 and 27-35% 
for HD 191612): this means that the presence of 
a high-temperature component is not the only 
reason for the overluminosity. 

Finally, we also analyze the X-ray lightcurves. 
First, one ought to note that the ROSAT data 
and Einstein count rates of HD 148937 agree 
well with the spectral properties derived from the 
XMM-Newton observations (Fig. [T2J). Second, 
lightcurves were derived from the XMM-JVewion 
datasets for a large range of energy domains and 
time bins. They were subsequently analyzed by 
X 2 and pov tests ( Sana et al. 2004 ) but no sig- 
nificant variation was found. It therefore appears 
that the flux of HD 148937 at X-ray energies is con- 
stant, within the error bars, on short-term ranges 
(over the duration of an observation) as well as on 
timescales of decades. 

5.1. High resolution (RGS) 

The high-resolution (RGS) X-ray spectrum of 
HD 148937 is presented in Fig. Q31 Although 
the signal-to-noise ratio is limited and in no ways 
comparable to that obtained for RGS spectra of 
closer objects, it is nevertheless possible to get 
some information from it. First, a global fit 
was undertaken: the lines were fitted by trian- 
gular profiles, which are particularly suited to get 
an idea of their width and their decentering (if 
any), while the residual continuum was fitted by 
a bremsstrahlung model. Only the amplitude of 
the triangles was allowed to vary from line to line, 
the overall shape being uniform. The best fit re- 
sults in a line center position of — 195±707kms~ 1 , 
a red edge at 1163 ± 512km s _1 , and a blue edge 
at -2327 ± 657kms _1 , resulting in a FWHM of 
1745 ±416kms _1 . The X-rays lines are thus quite 
large and do not present significant blue/redshift, 
but their profiles are slightly skewed. In this con- 
text, it might also be interesting to note that the 
Neix/Nex ratio f or this 06 s tar fit s nicely in the 
sequence found bv lWalbornl (|2007l) . 

To get more precise information, we decided to 
focus on the O Vin A 18.97 Lya line, the strongest 
feature in the spectrum that is rather free from 
blends with other lines. The line is broadened 
and has a FWHM around 2500kms _1 . A binned 
version of the data suggests that the line might 
display some structure (which could actually also 
be present on the Ly/3 line, see below). We have 



attempted to fit the Oviii Lya line with a ex- 
ospher ic line profile model following the formal- 



ism of Kramer et al. I (J2003J). The main assump- 
tions are that the X-ray emission originates from 
material distributed throughout a spherical wind, 
above a radius r > Rq > i?» and that the hot 
plasma follows the bulk motion of the cool wind. 
Doppler broadening due to macroscopic motion 
hence provides the main source of line broaden- 
ing. The line emissivity is assumed to scale as 
e oc p 2 r~ q where the r~ q term accounts for a ra- 
dial dependence of the filling factor of the X-ray 
plasma. The free parameters of this model are 
thus i?0j 9 (that we take equal to zero here) and 

Ta <* = 4 nv M R ' tr± e cnarac teristic optical depth 
at wavelength A. For the terminal wind velocity, 
we have first adopt e d = 2285kms _1 as derived 
bv lHowarth et al. I (|l997h . The best fit to the un- 
binned line profile is obtained for Rq ~ 2.25 i?* 
and ta,* = 0. However, larger values of ta,* are 
also possible for larger i?o- The fits to the binned 
profiles roughly confirm this picture, but bring up 
another (actually deeper) minimum in the \ 2 con- 
tours around Rq ~ 1.65 i?* and ta,* = 0.4. While 
the unbinned data hence favour a flat-topped pro- 
file, the binned spectra rather suggest a skewed 
profile. In any case, we note that the fits suggest 
a rather broad line. 

We have applied the same procedure to the 
O vin A 16.01 Ly/3 line. The best fit is obtained for 
Rq ~ 1.85 i?* and ta,* = 0, in reasonable agree- 
ment with the results obtained for the Lya line. 
Note that the strong Nex A 12.13 Lya line is prob- 
ably blended with another line on its red sidqj and 
can thus not be fitted easily with our model. 

We have repeated the fits of the O vm A 18.97 
Lya and A 16.01 Ly/3 lines with the same exo- 
spheric model and the same parameters as be- 
fore, except for the terminal velocity where we 
adopted Woo = 2600 km s -1 instead (see section 4.3 
above). The shape of the % 2 contours is essentially 
the same as found before with Voo — 2285 km s -1 , 
but, as can be expected from the increase of the 
wind velocity, the minimum is shifted to some- 
what lower values of Rq. In fact, the best fit 



3 There are a number of lines from various iron ions (e.g. 
Fe xxm A 12.193, Fexvn A 12.264 and Fexxi A 12.286) that 
could possibly be responsible for this blend. Given the 
presence of other strong lines of this ion, Fe XVII is probably 
the best candidate. 
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to the unbinned Lya line profile is now obtained 
for R ~ 2.05 i?* and r A ,* = 0. The fits to the 
binned profiles now yield the lowest % 2 around 
Ro ~ 1.45 R* and r A ,» = 0.0. For the Ly/3 line, 
the best fit remains at Rq ~ 1.85 i?* and t a ,* = 0. 

Over recent years, the He-like triplets, consist- 
ing of a forbidden (/), an intercombination (i) 
and a resonance (r) lines, of various ions have 
been used as plasma d iagnostics for a number 
of O-type stars (see e .g. Leutenegger et al.~ll200d 



Oskinova et al. 20061 ). In fact the H = f/i ra- 



tio has been shown to be a sensitive diagnostic 
of the dilution of the UV radiation field in the 
line emission region (e.g. lPorquet et al. llioOlh . In 
the case of HD 148937, the only He-like triplet 
with a reasonable level of exposure is the Neix 
triplet at AA 13.447 (r), 13.548 + 13.551 (i) and 
13.697 (/). However, even for this complex the 
data do not allow us to perform a quantitative fit 
of the line strengths. Still, the data show that 
the 1Z ~ 1 (see Fig. [14]). Hence this ratio must be 
below its co llision equilibrium low-de nsity limit of 



IZo = 3.1 ( Leutenegger et al. 20061 ). Note also 



that the / component could be blended with un- 
resolved Fexix lines at AA 13.73 - 13.74. These 
features could lead to an overestimate of the ac- 
tual strength of the / line. The reduced strength 
of the / component of th e Neix triplet is quit e 



typical for O-stars (see e.g. Oskinova et al. 2006h . 



6. Conclusions 

With its strong Cm A 4650 lines, HD 148937 is 
a true Of?p star. However, it is unclear if it con- 
stitutes a perfect 'twin' of HD 108 or HD 191612. 
This paper reports a first thorough, multiwave- 
length variability study of HD 148937. 

In the visible domain, its photometry does not 
change much, and its S-Dor status can not be con- 
firmed. The visible spectrum does vary, but the 
changes are limited to the main lines formed in 
the stellar wind, i.e. the Balmer Hydrogen lines 
and Hen A 4686. This variability occurs with very 
small amplitudes and an analysis of the changes in 
the Ha line reveals a periodicity of 7.031±0.003d. 
However, the He I and Cm lines seem constant. 
In addition, no large (>10kms _1 ) variations of 
the RVs could be identified and no other typical 
signature of binarity (e.g. blended spectrum) was 
detected. Finally, model atmosphere fits yield pa- 



rameters typical of an 05-6V/III star; they also 
suggest the presence of a non-spherical component 
to the stellar wind. 

In the X-rays, the spectrum appears thermal in 
nature, with a dominant cool component (0.2keV), 
broad unshifted X-ray lines, and an order-of- 
magnitude overluminosity compared to normal 
O-type stars. It looks nearly identical to those of 
HD 108 and HD 191612, though with a larger flux 
(because of the smaller distance) and thus a bet- 
ter signal-to-noise ratio. No significant short-term 
or long-term variation of the X-ray flux could be 
brought to light. 

Three stars in our Galaxy share a very pecu- 
liar spectral characteristic, the presence of strong 
C in lines. Is HD 148937 completely similar to the 
other two? Using only the X-ray results, the an- 
swer would be clearly yes. From the visible spec- 
troscopy, however, the answer is unclear. The 
variability of HD 148937 could either be related 
to the small-scale changes or to the large varia- 
tions seen in HD 108 and HD 191912. Consider- 
ing the former, additional data are needed to find 
if these small-scale changes are truly periodic for 
the other Galactic objects, and therefore deter- 
mine the origin of this short-term, small-scale vari- 
ability. In this case, it also remains to be seen if 
HD 148937 presents large variations of its H, Hei, 
and Cm lines - from our dataset, this can only 
occur on very long timescales (likely even larger 
than the 55yrs of HD108). On the other hand, 
the small-scale changes of HD 148937 could repre- 
sent a scaled version of the phenomenon observed 
in HD 108 and HD 191912, thereby explaining the 
lack of variations of the He I and C in lines. In this 
case, the 538d period of HD 191612 and the 55yrs 
timescale of HD 108 would here be replaced by a 
much shorter, ~7d period. 

To ascertain the origin of this variability and 
constrain more its properties, additional data are 
clearly needed: short-term monitoring of HD 108 
and HD 191612, long-term observations at high 
resolution and spectropolarimetry of HD 148937. 
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Fig. 1.— Hipparcos photometry of HD 148937. 



Fig. 2. — Visible spectrum of HD 148937 (median 
combination of our four FEROS spectra). The 
most important lines are labeled and tickmarks 
are drawn at their rest wavelengths. 

Fig. 3 — RVs and EWs of the HciA4471 and 
HellA4542 lines, and the ratio of the EWs. 
SMARTS data are presented by filled circles and 
for the right panel, FEROS by stars, UVES by a 
cross (shown at HJD+1000.), UCLES by an open 
triangle (shown at HJD+4000.). The limits for 



the spectral types come fro m IConti fc Alschuler 
(119711) : IConti fc Frostl (Il977l) . 



Fig. 4.— RVs and EWs of Ha and the DIB at 
A6613.62A, shown on the same scale (SMARTS 
data only). 



Fig. 10. — The H/3 line profile in the spectrum of 
HD 108 appeared particularly variable during our 
observing run of 2001: the thin solid line shows 
the spectrum taken on Sept. 10, the thick solid 
line the spectrum taken on Sept 11, dotted lines 
the spectra taken on Sept. 13 and 15, and dashed 
lines the spectra taken on Sept. 18. 



Fig. 11.— EPIC X-ray spectrum of HD 148937 
with the best-fit 3T model. The EPIC-pn data 
(drawn in green) appear at higher ordinates that 
those of the two EPIC-MOS (shown in black and 
red for MOS1 and 2, respectively). This figure 
only appears in colors in the electronic version of 
the journal. 



Fig. 5. — TVS and mean low-resolution spectrum 
of HD 148937. 

Fig. 6. — Left: long-term variability of the Ha line 
in the 2006 low-resolution spectra (from month to 
month) . Similar figures can be drawn for the other 
years of the campaign. Right: variability of the 
H/3 line in our FEROS spectra (2003: dotted line 
and black solid line for observations taken on May 
24 and 25, respectively; 2004: thin dashed line, 
2005: thin solid line). 

Fig. 7. — Top: Periodogram derived from the 
SMARTS spectra covering the Ha line (see text). 
Bottom: Spectral window of the observations. 



Fig. 12.— Normalized count rate: the ROSAT 
and Einstein count rates were divided by the 
count rate expected from the best-fit 2T model. 
A value of 1 indicates that no flux change has oc- 
curred. 



Fig. 13. — Binned (3 channels per bin), com- 
bined (RGS1+2, 1st and 2nd order) X-ray spec- 
trum of HD 148937. The most important emission 
lines are labeled. Note that the spectrum was not 
smoothed. 



Fig. 8. — RVs and EWs of Ha phased with a 
period of 7.031d (To was arbitrarily chosen as 
HJD=2 452 680). 

Fig. 9. — Comparison of the average FEROS spec- 
trum of some lines (solid black line) and their 
best fit CMFGEN model (dash-dot red line, for 
log(M) = —7.5 and the parameters listed in Sect. 
4.3). This figure only appears in colors in the elec- 
tronic version of the journal. 



Fig. 14.— The Ne IX fir He-like triplet in the RGS 
spectrum of HD 148937. The rest wavelengths 
of the different components are indicated by the 
labels. The / component is clearly suppressed 
compared to the collision equilibrium low-density 
value of IZq — 3.1. 
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Table 1: Mean and dispersion of the measured RVs and EWs. 



Line 


EW interval 


KV smarts 


EW smarts 


KV/ eros 








(A) 


(kms -1 ) 


(A) 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 


He ii A 4199.83 


4193.83-4205.83 


-46.0±12.2 


0.495±0.054 


-38.0±1.4 


-40.3 


-41.3 


He I A 4471.512 


4465.512-4475. 


-59.3±10.4 


0.325±0.034 


-41.0±8.1 


-30.7 


-51.3 


He II A 4541.59 


4535.-4548.18 


-38.7±7.5 


0.683±0.035 


-32.2±2.0 


-32.0 


-33.1 


N in A 4634.25 


4631.-4644.25 


-24.4±7.3 


-0.963±0.067 


-29.7±4.6 


-31.3 


-26.7 


N in A 4641.02 




-37.4±6.7 




-43.6±0.3 


-50.2 


-44.3 


Cm A 4650 


4644.25-4653.40 




-0.470±0.037 








He ii A 4685.682 


4679.682-4691.682 


-8.8±7.5 


-0.761±0.079 


-18.3±4.7 


-22.8 


-19.4 


H/3 A 4861.33 


4645.33-4877.33 


-191.2±11.9 


1.218±0.120 


-182.0±4.6 


-185.0 


-187.7 


Cm A 5695.92 


5689.92-5701.92 


-34.7±11.9 


-0.205±0.041 


-22.1±0.4 


-28.9 


-23.2 


DIB A 5780.45 


P* ^7^7 A AC r^o/ 1 A r~ 

5774.45-5786.45 


— 19.8±11.0 


0.515±0.045 


— 13.6±0.5 


— 15.1 


gap 


CivA5811.98 


5805.98-5817.98 


-49.3±13.6 


0.196±0.035 


-33.3±3.6 


-29.6 


gap 


Civ A 5812 (cor.) 




-42.7±8.6 










Ha A 6562.85 


6550.85-6574.85 


3.0±12.6 


-1.176±0.237 


-7.1±7.6 


-2.3 


-12.4 


Ha(cor.) 




2.2±12.9 










DIB A 6613.42 


6609.42-6617.42 


-2.2±10.4 


0.106±0.020 


-3.4±1.1 


-8.5 


-5.2 



Table 2: Best-fitting models and X-ray fluxes at Earth. 



Type 


10 22 cm" 2 


kTi 
keV 


normi 
lO^cm" 5 


JV 2 H 
10 22 cm" 2 


kT 2 
keV 


norm,2 
lO^cm" 5 


Nf 
10 22 cm" 2 


kT 3 
keV 


norm3 
10- 3 cm- 5 


x2(dof) 


f abs 

(in 10- 


funabs 

Jx 

12 ergcm~ 2 s _1 ) 


2T 


0.49° : |i 


n oq0.23 
u - zo 0.22 


33.9 3 |:* 


1 03 112 


1 34 1 - 37 


4.42H1 








1.34 (1804) 


3.0 


7.0 


3T cool 


n 0.06 


o.o88;81 


i9.5?l:? 


^0.57 


n 23°- 24 


38.4ii 


1 90I.33 


1-371:18 


4.52t: 3 l 


1.23 (1801) 


3.1 


8.6 


3T hot 


o.oi 


o.258 : 2 t 


i.57i;S 


35 39 


632 




53 ' 63 
u.oo 44 


2 -oo 2 :^ 


2.25 2 : 3 ? 


1.13 (1801) 


3.2 


7.8 
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